EPJ manuscript No. 

(will be inserted by the editor) 



Gravitational collapse of a magnetized fermion gas with finite 
temperature 



o 



00 

o 

cn 
p 

o 

cn 



I. Delgado Caspar^, A. Perez Martinez^, Roberto A. Sussman"^, and A. Ulacia Rey^'^ 

^ Institute de Geofi'sica y Astronomia (IGA), Calle 212 No 2906, La Lisa, La Habana, 11600, Cuba, e-mail: idelgado@iga. cu 
^ Institute de Cibernetica, Matematica y Fisica (ICIMAF), Calle E esq 15 No. 309 Vedado, La Habana, 10400, Cuba, e-mail: 
auroraSicmf . inf . cu 

^ Institute de Ciencias Nucleates, Universidad Nacional Autonoma de Mexico (ICN-UNAM). A. P. 70-543, 04510 Mexico D. 

F., e-mail: sussman@nucleares.unain.mx 
* Institute de Cibernetica, Matematica y Fi'sica (ICIMAF), Calle E esq 15 No. 309 Vedade, La Habana, 10400, Cuba, e-mail: 

alainSicmf . inf . cu 

Received: date / Revised version: date 

Abstract. We examine the dynamics of a self-gravitating magnetized electron gas at finite temperature 
near the collapsing singularity of a Bianchi-I spacetime. Considering a general and appropriate and phys- 
ically motivated initial conditions, we transform Einstein-Maxwell field equations into a complete and 
self-consistent dynamical system amenable for numerical work. The resulting numerical solutions reveal 
the gas collapsing into both, isotropic ("point-like") and anisotropic ("cigar-like") singularities, depending 
on the initial intensity of the magnetic field. We provide a thorough study of the near collapse behavior 
and interplay of all relevant state and kinematic variables: temperature, expansion scalar, shear scalar, 
magnetic field, magnetization and energy density. A significant qualitative difference in the behavior of the 
gas emerges in the temperature range T ~ lO'^K and T ^ lO^K. 

PACS. 4.20.-q, 04.40.4-b, 05.30.Fk 



1 Introduction 

Astrophysical systems provide an ideal scenario to exam- 
ine the effects of strong magnetic fields associated with 
self-gravitating sources under critical conditions. In such 
conditions, we expect non-trivial coupling between grav- 
itation and other fundamental interactions (strong, weak 
and electromagnetic), and from this interplay important 
clues of their unification could emerge. 

The presence and effects of strong magnetic fields in 
compact objects (neutron, hybrids and quark stars) have 
been studied in the literature (see [UISIISIIIIEIEIIZIIHI , and 
references quoted therein) , assuming various types of equa- 
tions of state (EOS) has been obtained and considering in 
some of these papers numerical solutions of the equilib- 
rium Tolman-Oppenheimer-Volkov (TOV) equation. 

As proven in previous work pilSll^ lTUlfTT] , the presence 
of a magnetic field is incompatible with spherical symme- 
try and necessarily introduces anisotropic pressures. How- 
ever, by assuming that anisotropics and deviations from 
spherical symmetry remain small, several authors [SlfTUl 
I12| have managed to compute observable quantities of ide- 
alized static and spherical compact objects under the pres- 
ence of strong magnetic fields by means of TOV equations 
that incorporate these anisotropic pressures. Moreover, it 
is evident that much less idealized models would result 



by considering the magnetic field and its associated pres- 
sure anisotropics in the context of TOV equations under 
axial symmetric (or at least cylindrically symmetric) ge- 
ometries [5]. 

As an alternative (though still idealized) approach, and 
bearing in mind the relation between magnetic fields and 
pressure anisotropy, we have examined the dynamics of 
magnetized self-gravitating Fermi gases as sources of a 
Bianchi I space-time [TSEHIIS] , as this is the simplest 
non-stationary geometry that is fully compatible with the 
anisotropy produced by magnetic field source. 

Evidently, a Bianchi I model is a completely inade- 
quate metric for any sort of a compact object, as all ge- 
ometric and physical variables depend only on time (and 
thus it cannot incorporate any coupling of gravity with 
spatial gradients of these variables). However, the use of 
this idealized geometry could still be useful to examine 
qualitative features of the local behavior of the magnetized 
gas under special and approximated conditions. Specifi- 
cally, we aim at providing qualitative results that could 
yield a better understanding of the conditions approxi- 
mately prevailing near the center and the rotation axis of 
less idealized configurations, where the angular momen- 
tum of the vorticity and the spatial gradients of the 4- 
acceleration and other key variables play a minor dynam- 
ical role. 
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The main objective of the present paper is to study 
the dynamical evolution of a self-gravitating magnetized 
Fermi gas at finite temperature, under the extremely crit- 
ical conditions near a collapsing singularity. In particular, 
we examine the possibility that a finite temperature may 
produce the dynamical effect of "slowing down" or revers- 
ing the collapse of the magnetized gas. Also, we aim at 
verifying how a finite temperature affects the evolution of 
other thermodynamic variables, such as the magnetization 
and energy density. This work constitutes a continuation 
and enhancement of previous work [13 . 

The paper is organized as follow. In section II we derive 
the EOS for a dense magnetized electron gas at finite tem- 
perature. In section III we lay out the dynamical equations 
for the evolution of our model by writing up the Einstein- 
Maxwell system of equation for the specific source under 
consideration. In section IV we introduce the local kine- 
matic variables, which are later used in section V to write 
the Einstein-Maxwell equations as a system of non-linear 
autonomous differential equations, rewriting it afterwards 
in terms of physically motivated dimensionless variables. 
The numerical analysis of the collapsing solutions and the 
discussion of the physical results are given in the section 
VI. Our conclusions are presented in sections VII. 

2 Magnetized Fermi gas as a source of a 
Bianchi I background geometry 

Homogeneous but anisotropic Bianchi I models are de- 
scribed by the Kasner metricQ 



field. Notice that the anisotropy in is produced by 
the magnetic field B. If this field vanishes, the energy- 
momentum tensor reduces to that of a perfect fluid with 
isotropic pressure (an ideal gas of electrons complying 
with Fermi-Dirac statistics). 

The EOS for this magnetized electron gas can be given 
in the following form [41 ITT1I16) : 

P\\ = -Q,-BNU + —, (5) 

Pi- = - ^, (6) 

U = -Qe + TS + fiN+—. (7) 

where Me = -- {dQe/dB) is the electron-system magneti- 
zation, S — — [dfle/dT) is the entropy, N = — {dfle/dix) 
is the particle number density and the term B^ /2 is Maxwell's 
classical magnetic contribution. We note that equations 
dS]) and ^ can be re-arranged as [HIH]: 

Pj_ = P\\ - BM, (8) 

The thermodynamic potential Q,, can be written as 
the summ: 

- QQ^'^ + (9) 

where we identify the vacuum term i?^^^, independent 
of the temperature and chemical potential, and the sta- 
tistical term, flf^^'^ . As is known, ^^^^ has non- field- 
dependent ultraviolet divergencies, hence after renormal- 
izing the Schwinger expression we obtain 117J : 



ds^ = ~dt^ + Qi [tf dx^ + Q2 {tf dy^ + Q3 {tf dz^, (1) 

so that spatial curvature vanishes and all quantities de- 
pend only on time. Assuming a comoving frame with coor- 
dinates x°' = [t, X, y, z] and 4- velocity tt° — 51 1 the energy- 
momentum tensor for a self-gravitating magnetized gas of 
free electrons is given by: 



Tfc'' = (C/ + P)u°Ub + PV + 7T,", P=p- 



2BM 



(2) 



where B is the magnetic field (pointing in the z direction), 
U is the energy density (including the rest energy of the 
electrons), M = - [dQ/dB) with = i7e + and Qe 
is the thermodynamical potential of the electron gas, P is 
the isotropic pressure and 7T(j° is the traceless anisotropic 
pressure tensor: 



77fc° = diag [0, iT, 77, -277] , n 



BM 



(3) 



We can write the energy-momentum Q tensor as: 

r,'^ = diag[-C/,Fx,P±,P||] , (4) 

which respectively identifies Pl and P\\ as the pressure 
components perpendicular and parallel to the magnetic 



-TO^s) X 



[esBf 



(10) 



^ Unless specified otherwise, we use natural units. 



X esB coth {esB) - 1 

V " / 

r2fQ^^(7?,T,^) = 

J-00 

where d{l) = 2— 610 is the spin degeneracy of Landau levels 
with / ^ 0, nie is the electron mass, ei = + 2\eB\l + m\ 
and /3 = 1/T. 

Finally, we remark that in a high density regime (which 
we shall consider henceforth) the contributions to the EOS 
of the vacuum term and the Maxwell term {i.e. f^^^^ and 
73^/2) can be neglected whenever the magnetic field is less 
than critical {i.e. B < Be — m^/e). Thus, in subsequent 
calculations the leading contribution to the EOS will come 
from flf^^^ (see Appendix A). 



3 Einstein-Maxwell equations 

Since we are interested in the critical relativistic regimen, 
the dynamics of the magnetized gas whose EOS we have 
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described in the previous section must be studied through the kinematics of local fluid elements through the covari- 
the Einstein field equations in the framework of General ant derivatives of it". For a Kasner metric in the comov- 
Relativity: ing frame endovifed with a normal geodesic 4-velocity, the 

only non-vanishing kinematic parameters are the expan- 
-Rgfii^ — nT^i, (12) ^^^'^ scalar, 0, and the shear tensor ctq 



Gn 



Rn 



together with the balance equations of the energy-momentum 
tensor and Maxwell's equations, 







(Jap 



-F^" = 0, 



= 0, 
= 0, 



(13) 
(14) 



where k = SttGn and Gat is Newton's gravitational con- 
stant, while square brackets denote anti-symmetrization 
in pi)) . 

Assuming absence of annihilation/creation processes, 
so that particle numbers are conserved, leads to the fol- 
lowing conservation equation: 



0, 



Nu", 



(15) 



where es the particle number density. From the field 
equations (fT^ we obtain: 



Gx 



Q2Q3 
Q2Q3 
y _ Q1Q3 
^ ~ Q1Q3 
Q1Q2 



Q1Q2 



Q2 
Qi 
Qi 
Qi 
Qi 
Q1Q3 



Q3 
Q3 



W2 
Q2 
Q2Q3 



Q1Q2 QiQz Q2Q3 



kU. 



(16) 
(17) 
(18) 
(19) 



(23) 



where h^p = UaUp + gap is the projection tensor and 
rounded brackets denote symmetrization. These parame- 
ters take the form: 



Q _ Qi ^ Q2 ^ Q3 
Qi Q2 Q3 



(24) 



a^" = diag [a/ , a/ , a/ , O] = diag [^1,^2,^3,0], (25) 
where: 



2Qa 
3Qa 



IQt 



IQe 

" 3Qc' 

a ^ b ^ c (a, 6, c 



1,2,3). (26) 



The geometric interpretation of these parameters is straight- 
forward: represents the isotropic rate of change of the 
3- volume of a fluid element, while describes its rate of 
local deformation along different spatial directions given 
by its eigenvectors. Since the shear tensor is traceless: 
= 0, it is always possible to eliminate any one of the 
three quantities {Si, S^) in terms of the other two. 
We choose to eliminate Si as a function of {S2, S-^). 



where Q = Q-au" — Q^t- From the conservation of the 
energy-momentum tensor (|13p we obtain: 



5 Dynamical equations 



Maxwell's equations ((T4)) yield 



Qi Q2 2B 



(21) 



and from the particle number conservation (|15l) leads to: 

(22) 



4 Local kinematic variables 

Einstein-Maxwell field equations are second order system 
of ordinary differential equations (ODE's). In order to 
work with a first order system of ODE's, it is useful and 
convenient to rewrite these equations in terms of covariant 
kinematic variables that convey the geometric effects on 



By using equations (1^^ and (|26l) we can re-write the sec- 
ond derivatives of the metric functions in p^ . (fTTl) y ([T^ 
as first order derivatives of 0, S2 and S3. After some alge- 
braic manipulations it is possible to transform equations 
((T6| - ((22|) as a first order system of autonomous ODE's: 



S2 
S3 



3 ' 
2k 



P^-P^^)-0S2, 



0S. 



B = 2B [ S-i- -0 



f 



Dct 

1 

Det 



[hN,T ~ f2U,T] , 



together with the following constraint: 



kU 



(S2 



SyS-t 



02 

3 ' 



(27) 
(28) 
(29) 
(30) 
(31) 

(32) 
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where: 



Det = C/,^7V,T - U,tN,^, 
h - (^3 - ^) (P± + U) 



e 



(33) 



^3 + y ) (^[| +U)- U^bB, (34) 



f2 = 0N + N,bB. 



(35) 



We introduce below the following dimensionless evolu- 
tion parameter, 

d _ I d 



Be rUe 

and write the EOS in the form: 

u = x^u{p,^i,<p), 

N - (A/me) r„ if], fi, 



T 
rrie. 



(37) 



(38) 
(39) 
(40) 
(41) 



where A = STrm,,/ (X^) with Ac the Compton wavelength 
of the electron and Hq is a constant with inverse length 
units that sets the characteristic length scale of the sys- 
tem, which we have chosen as 3-ffQ = kA Hq = 0.86 x 
10-1° m-i, so that I/Hq = 1.15 x 10^° m is of the or- 
der of magnitude of an astronomic unit. It indicates that 
our simplified model is examined on local scales smaller 
than cosmic scales. In cosmological sources and models 
[TF iJo = 0.59 X 10-26 m-i would play the role of the 
Hubble scale constant, this value is a much greater length 
scale. The functions 5*2 and S'3 are the components of the 
shear tensor normalized with this scale, while r is the di- 
mensionless time. 

Substituting (I36l)-(I1II) into the system ([27l)-(l3ll) we 
obtain: 



5*2, r — r± — r'li — 3HS'2, 
b^r = 26(53-27^), 

1 



Det 



while the constraint (15^ becomes 



SFif — ~<S'2 — So — <S'2>5'3 + 37i^, 



(42) 
(43) 
(44) 

(45) 
(46) 



(47) 



and the auxiliary parameters of the previous system take 
the form: 

/i = (^3 - 2n) (rx - 2ru,bb) - {S3 + n) ^ - srun, 

/2 = mrN + 2rjv,f, {S3 - 2H) , (48) 
Det — ru^j2rN,4, — ru^^rN^f, . 



These differential equations form a complete and self- 
consistent system whose numeric integration fully deter- 
mines the variables H, 5*2, 5*3, 6, p,, (p, thus allowing us to 
study the dynamical evolution of a local volume element 
of a gas of magnetized electrons with finite temperature. 
In particular, all thermodynamical magnitudes, such as 
the particle number density, the magnetization, the energy 
density, the entropy and the pressures, are functions of 
^(t), //(t), (/)(t) and thus can be determined from the nu- 
merical solution of this system. On the other hand, H, S2, 
S3 provide all the necessary information on the kinematic 
evolution of volume element, and in particular, its proper 
volume and the metric functions. 

Since V = — det g^ff = Q1Q2Q3, we obtain by means 
of and the local volume in terms of H: 



V{t)= V (0) exp 3 / HdT 



(49) 



T = 



where we remark that the sign of H (r) implies expansion 
if H (r) > 0, and collapse if TL (r) < 0. Besides this point, 
equations (|24p and (1^51) lead to: 



{H + S,) dT 



Qi{t) Qi(0)exp 



T = 



(i = 1,2,3). 
(50) 

where Si = — (^2 -|- ^3). In the following section we under- 
take the numerical study of the system (|^ - (|Tf)) . focusing 
specifically on the the collapsing regime. 



6 Numeric analysis and physical 
interpretation 

For the numerical study we assume a magnetized elec- 
tron gas at high density: /ie(0) = 2, which means that 
the chemical potential is 2me. The initial values of the 
magnetic field and temperature were chosen in the ranges 



6(0) - 10-5 to &(0) - 10-4 and 



10-" to 



10-^ 



respectively 0. Together with "H (0) < 0, we consider the 
conditions 5*2(0) — 0, ^3(0) = 0, -1-1, which correspond to 
the cases with zero initial deformation and initial defor- 
mation (shear) in the direction of z axis, respectively. The 
calculation has been done using the fourth-order Runge- 
Kutta method with the local truncation relative error less 
than 10-6. 

The numerical solutions for the function % for the as- 
sumed values 'H(O) show that % — —00, which implies 
that the volume element evolves to a singularity (see equa- 
tion (j49p ) . This is exemplified in figure [TJ where numeri- 
cal solutions are displayed for the expansion scalar. These 
curves correspond to different values of the initial tem- 
perature in the range (/)(0) = 10-^ to (/)(0) = 10-^, fixing 
the rest of the initial conditions on the values /i(0) = 2, 



^ For instance B ~ lO^G (b ~ IQ-^), T ~ lO^K {(j) ~ IQ-^) 
and ^ ~ IMeV {jl — 2) correspond to the physical parameters 
of a white dwarf type astrophysical object. 
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Fig. 1. Numerical solutions for the expansion scalar H (r). 
The box in the lower corner is an amplification of the graphic. 



6(0) 5 X 10-^ S2{0) = and 53(0) = 1. Notice that 
in this regime (as given by these initial conditions) where 
the high densities are dominant, we do not obtain a di- 
rect relation between the values of initial temperatures 
and the collapse time. However, the collapsing time di- 
minishes when the values of the initial magnetic field in- 
creases, which agrees with the results obtained in [T3], 
where T = was assumed. 

In all examined configurations the strength of the mag- 
netic field diverges during the collapse of the volume el- 
ement, an expected development due to compression of 
the magnetic fluid. Since the energy-momentum tensor 
(m takes the perfect fluid form (P|| — P±) for zero mag- 
netic field, we can identify the magnetic field as the factor 
introducing anisotropy in the dynamical behavior of the 
fermionic gas. In particular, this anisotropy in the pres- 
sures (P[[ 7^ P±) must yield different evolution in different 
directions, which must be evident in a critical stage such as 
the collapsing regime. When the evolution occurs in such 
a way that all metric coefficients tend to zero, the singu- 
larity is necessarily isotropic "point-like" , whereas when 
two metric coefficients evolve to zero and the third one 
to -l-cxD the result is an anisotropic singularity, "cigar-like" 
(see the definitions of these types of singularities in [H]). 
It follows from (jSOl) that the evolution of the terms Si + H 
to ±00 implies that the metric coefficient Qi evolves to 
cither +00 or 0, which determines the type of singularity. 
Intuitively, we expect an isotropic point-like singularity if 
the pressure is isotropic, as pressure diverges in all direc- 
tions, but a large pressure anisotropy (which necessarily 
corresponds to large magnetic field) should lead to a qual- 
itatively different direction dependent critical behavior of 
the pressure, which should result in an anisotropic cigar- 
like singularity characterized by the divergence of only the 
pressure parallel to the magnetic field. 



Fig. 2. The plots of the functions Si -fH for i = 1, 2, 3. The 
system has an initial deformation (shear) in the same direction 
of the magnetic field: z axis. The collapse is in the form of a 
"cigar-like" singularity in the 2 direction. 



S^{0)=0 S3 (0)=1 




' 1 ' 




- - - (0)=S3 (0)=0 








p = p 

1 II 










4.96x10* 


] 






4.95x10* 


✓/ 1 




4.95x10* 4.96x10* ! 

1 




— * -*- 











1x10* 2x10* 3x10* 4x10* 5x10* 

Fig. 3. Pressure parallel and perpendicular to the magnetic 
field. We plot P||/A vs P^/X for initial shear 5*2(0) = 5*3(0) = 
and 52(0) = QA53(0) = 1 taking (0) = 10"^ A = 8Tvme/{\i). 
The case 6 = is shown for reference. 



In the numerical calculations whose initial conditions 
assume zero shear: 5*1(0) — 82(0) = 6*3(0) = 0, we al- 
ways obtain a point-like singularity, independently of the 
selected values for the remaining initial conditions. How- 
ever, when the initial deformation (shear) is positive in 
the direction of the magnetic field {z axis), we obtain a 
cigar-like singularity. Figure [5] displays the solutions of 
S + nior initial conditions: /i(0) = 2, 5(0) = 5 x 10-^ 
5*3 (0) = 1, 5*2 (0) = and (/)(0) = lO^^, 10-^, lO^^. The 
nine curves displayed in the figure correspond to three sets 
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of functions: Si + H [i — 1,2,3), each one correspond- 
ing to three different values of the temperature: </> (0) — 
10~^, 10~^, 10~^. As we mentioned previously, in these 
cases the gas collapses into a cigar-like singularity in the 
direction of the magnetic field, and when the initial shear 
is zero the three curves Si+H tend to — oo for each initial 
value of temperature. 

In figure [3] we compare the behavior of the anisotropic 
pressure in the directions parallel vs perpendicular during 
the collapse. The dashed line corresponds to the case with 
zero initial shear, while the star-gray line corresponds to 
the case of initial positive shear in the z direction. In 
both cases the remaining initial condition were fixed as: 
0(0) = 10-3, fi{0) = 2 and 6(0) = 5 x 10"^ As shown 
by the figure, the pressures exhibit the expected behav- 
ior: for initial zero shear we have P|| ~ Pi and a point 
singularity, but for initial positive shear in the z direction 
we have a highly anisotropic evolution P|| ^ P±, which 
signals a cigar-like singularity. The case P|| = P± with 
zero magnetic field is shown for comparison (solid line). 

For initial values in the range 0(0) ~ 10"^ — 10~^ the 
temperature decreases as the collapse proceeds, while for 
0(0) ^ 10^"^ and higher values, it increases. This behavior 
of temperature can be explained by considering that for 
small initial values {4>{0) ^ 10~^ — 10~^) the effects of the 
magnetic field (alignment of electron spins with the mag- 
netic field) predominate over the effects of temperature 
(random motion) throughout the whole evolution. There- 
fore, we argue that the system "cools down" (tempera- 
ture decreases) during the collapse process as the mag- 
netic field increases. On the other hand, for high initial 
temperatures (0 (0) = 10""^), the system is "heated" (tem- 
perature increases) during the collapse, which shows the 
predominance of the thermal effect of temperature over 
the magnetic field effect. 

We display in figure 0] two dashed curves for the evolu- 
tion of the temperature from the initial value (0) — 10^^, 
together with solid curves that depict temperatures start- 
ing at 0(0) = 10^^. In both cases the initial conditions 
correspond to zero initial shear and positive deformation 
in the z axis. The initial values of magnetic field and chem- 
ical potential were fixed at 6(0) = 5 x 10^^ and //(O) = 2 
respectively. Notice that the symmetric configuration de- 
lays the collapse (yields a longer collapsing time), which 
may indicate a connection with stability of local fluid el- 
ements in compact objects. This could provide an impor- 
tant clue on the stability of a compact object made of a 
dense magnetized gas. However, verifying this possibility 
is beyond the scope of this article. 

We examine now the behavior of the magnetization 
and the energy density of the system during the collapse. 
We remark that the r-depending functions: 6, fl, fol- 
low from the numerical solution of the system (|42l) - (|46|) 
and (I47p . which allows us to compute the magnetization 
M (6, p,, 0) and U (6, p,, 0) for all values of r. The figures 
figures [5] and [S] have been done with the initial conditions: 
6 (0) = 5 X 10-^ S2 (0) = S3 (0) = 0. 

Figure [5] depicts the behavior of the magnetization vs 
the magnetic field. Notice how the gas exhibits the de 
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Fig. 4. Temperature of the system. The left hand side of the 
graph displays the temperature values (0i) corresponding to 
the two solid curves. The right hand side displays the temper- 
ature values (02) corresponding to the two dashed curves. 




Fig. 5. Magnetization versus magnetic field. The dashed curve 
(the solid curve) represents the magnetization taking (0) — 
10"^ (0(0) = 10-"^). 



Haas- van Alphen oscillations for small values of the mag- 
netic field and temperature. These oscillations have been 
observed in many magnetized systems, and are a direct 
consequence of the Pauli exclusion principle and the dis- 
creteness of the spectrum [50]. When the gas starts from 
initial values (0) = 10~^ the magnetization has an oscil- 
latory behavior, with these oscillations disappearing as the 
magnetic field increases. As a contrast, for initial values 
0(0) = lO^'^, the system is heated as the collapse pro- 
ceeds and the magnetization has a non-oscillatory behav- 
ior. These results, which were obtained in a time-dependent 
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*1 

0.0 0.4 0.8 1.2 




3.0x10"' 6.0x10' 9.0x10"' 1.2x10"' 



*2 

Fig. 6. Energy density versus temperature. The temperature 
values of the curve draws with solid line should be observed in 
the lower axis {(j>2) and in the top axis {(j>i), those of the curve 
drawn with dashed line. 



system, agree "with those obtained in pTll^ , -where the au- 
thors also studied the magnetization of an electron gas at 
finite temperature. This agreement sho^ws how the time 
evolution that we have studied preserves the properties of 
the same system when it is time-independent. 

In all the studied cases the energy density increases 
as the collapse proceeds, though this fact does not im- 
ply that temperature increases. As we have seen above for 
(?!>(0) ~ 10^^, the system is heated during the collapse, 
which yields a monotonic behavior of the energy density 
with respect to the temperature. On the other hand, for 
(j) (0) = 10^^ — 10"^ the system evolves in such way that 
temperature decreases in time as energy density increases, 
without the observed monotonic behavior of the energy 
density vs the temperature. In this latter case the chemi- 
cal potential and the particle density increase, even if the 
temperature decreases and the magnetic field increases, 
Hence, the energy density also grows during the collapse. 
Figure [6] displays two curves of the energy density vs tem- 
perature, with the dashed and solid line respectively cor- 
responding to 0(0) = 10-^ and 0(0) = 10"^. 

It is important to remark that the study we have un- 
dertaken here is based on an equation of state and a ther- 
modynamical potential {f2e) that come from a one- loop 
approximation. This approximation may not be appropri- 
ate in regimes in which the magnetic field overcomes its 
critical value. This issue deserves a separate study, pos- 
sibly in the context of non-perturbative calculations at 
high magnetic fields, which would mean a different ap- 
proach that may be applicable to the magnetized gas that 
we have examined here (see [53] and references quoted 
therein) . 



7 Conclusion 

We have examined the dynamical and thermodynamical 
behavior of a magnetized, self-gravitating electron gas 
at finite temperature, taken as the source of a simpli- 
fied Bianchi I space-time represented by a Kasner met- 
ric, which is the simplest geometry that allows us identify 
the magnetic field as the main source of anisotropy. We 
regard this configuration as a toy model that roughly ap- 
proximates a grand canonical subsystem of a magnetized 
electron source in the conditions prevailing near of the cen- 
ter and rotation axis of a compact object (in which spa- 
tial gradients of physical and kinematic variables may be 
regarded as negligible). The resulting Einstein-Maxwell 
field equations were transformed into a system of non- 
linear autonomous evolution equations, which were solved 
numerically in the collapsing regime for a chemical poten- 
tial: p, = 2me, a magnetic field in the range B ~ 10^ — lO^G 
and temperatures T ~ IC* — lO^K. 

For all initial conditions that we considered the gas 
evolves into a collapsing singularity, which can be (de- 
pending on the initial conditions) isotropic ( "point-like" ) 
or anisotropic ("cigar-like"). In aU the studied cases the 
isotropic configuration is the most stable, as the collaps- 
ing time is larger. This result may be connected with the 
stability of volume elements in less idealized configura- 
tions, an issue that is outside the scope of this paper and 
deserves a proper examination elsewhere. 

The collapse time decreases as the initial magnetic field 
increases, but we did not find a direct proportionality re- 
lation between this time and the initial temperature (this 
may be a consequence of having assumed a high density 
regime) . 

The behavior of the temperature as the collapse pro- 
ceeds also depends on initial conditions: for initial tem- 
peratures T ~ lO'' — lO^K the temperature decreases, 
while for values T > lO^K the temperature increases. 
As mentioned previously, this difference in the temper- 
ature behavior can be explained by the predominance of 
the magnetic effects (alignment of electron spins with the 
magnetic field) over the effects of temperature (random 
motion) during the evolution of the gas in the regime char- 
acterized by low temperature and strong magnetic field. 

The behavior of the magnetization and energy den- 
sity was studied by means of the numerical solutions of 
the system. We found a monotonic relation between the 
magnetization and the magnetic field for high tempera- 
ture values, but this relation does not occur for low tem- 
perature values. For small temperature values, when the 
magnetic field is also small, the electron gas exhibits the 
Haas- van Alphen oscillations. These results, obtained in a 
time-dependent system, agree with those obtained in j21[ 
122] for a time-independent system. 

In all the studied cases the energy density increases 
as the collapse process, even when the system is cooled. 
This occurs because the particle number density and the 
chemical potential increase during the collapse, causing 
the energy density to grow even when the temperature 
drops. 



8 



I. Delgado Caspar et al.: Gravitational collapse of a magnetized fermion gas with finite temperature 



The study we have presented can be readily applied to 
examine hadronic systems (complying vifith suitable bal- 
ance conditions and adequate chemical potentials). The 
methodology we have used can also serve as starting point 
to study the origin and the dynamics of primordial cosmo- 
logical magnetic fields. These potential extensions of the 
present work are already under consideration for future 
articles. 
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A Contributions to EOS 

The thermodynamical potential of an electron gas has two 
contributions, see equation Q, and thus we can express 
the magnetization as = Msqft + Mqpt- Taking it 
into account we can rewrite the EOS (P^ - (|Tfl) in the fol- 
lowing form: 

P|l = -{2^QPT _ B Msqft -^"^^^ ~ BMqft + — 



^P\\[B,fi,T) + AP\\{B), (51) 

P± = —^SQFT —^QFT TT 



= P^[B,ii,T) + AP^{B), (52) 

B^ 

U = -QSQFT + TS + ^iN -Qqft + -Y 



= U{B,^i,T) + AU{B), (53) 

where we identify by P±, P\\ and U the terms correspond- 
ing to our approximate EOS, while ^^"11, AP± and AU 
represent those that have been neglected. Note that these 
last terms are only field-dependent. 

In order to analyze the contribution of AP\\ , AP± and 
AU to its respective EOS, we evaluate the expressions 
P||/Z\P||, T±/AP± and U/AU in the sets of solutions ob- 
tained previously. In all the cases the numerical results 
show that these terms are not significant for values of 
magnetic field less than or equal to the critical magnetic 
field. 

To exemplify that mentioned above, in the figure [7] 
we show the fractions of the EOS neglected during the 
evolution of the gas starting from (f){0) = 10~^, b{0) = 
5 X 10-^ n{0) = 2, 82(0) = and ^'3(0) = 1. As we can 
observe these values remain near to zero for B < Be- The 
graph of magnetic field was included for comparison. 



Fig. 7. In the upper (lower) part are shown the fraction of 
contribution to EOS of the neglected terms (values of magnetic 
field) versus proper time. 
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